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A new ab initio theoretical formulation to calculate Z e ff and hence the positron annihilation 
rates is presented using the on-shell and half-offshell T-matrix scattering amplitudes without any 
explicit use of the scattering wave function. The formulation is exact and universal, as long as the 
dynamics of the systems can be described by a Lippmann-Schwinger type equation. It could serve 
as an effective tool as all the T— , K-, and S— matrix formulations, yield directly the scattering 
amplitudes; not the wave function. We also present a basic numerical test of the formulation. 
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I. INTRODUCTION 



Positron, being an antiparticle, interacts intimately (short-range interaction) with atomic electrons due to the 
absence of any restriction imposed by the Pauli Exclusion principle. Thus their annihilation studies, namely, the 
Positron annihilation spectroscopy and Positronium annihilation lifetime spectroscopy have emerged as two front line 
research areas, as they are expected to provide a more detail account of their close interactions with the target and 
bear the potential of various modern technological applications [1-3]. 

Up until now, the theoretical study of annihilations rate requires the evaluation of the scattering wave function. 
However, it is of fundamental importance to note that most scattering calculations (T-, K-, S-matrix) yield directly the 
scattering amplitudes; not the wave function. Thus, the evaluation of Z e f / normally requires a separate calculation for 
the wave function using these scattering amplitudes (like those in ref. [4]), or need to adopt of a different methodology 
that will yield the scattering wave function directly. Ryzhikh and Mitroy [4] have used the T-matrix formulation to 
evaluate the scattering wave function which has subsequently been used to evaluate Z e ff. Gribakin [5] has used an 
approximate form of the T-matrix equation to evaluate Z e ff. Here, we present a new ab-initio theoretical formulation 
based on the integral equation formalism, whereby Z e ff and hence the annihilation rates can be exactly calculated 
using on-shcll and half-offshell T-matrix scattering amplitudes but without any explicit use of the scattering wave 
function. 

Theoretically, the positron annihilation rates (A) are expressed in terms of Z e ff, the effective number of target 
electrons available to the incoming positron as [6-8] : 

A = nr^cZeffN sec~l (1.1) 

And Z e ff is defined in terms of the scattering wave function \ip£) as [6-8] 

JV 

Z e// (k) = (V+|^^(r j -x)|V+) (1.2) 

3=1 

where ro is the classical electron radius; c is the speed of light; Airr^c is the non-relativistic electron-positron annihi- 
lation rate; irr^c is the same for the spin-averaged case of two-gamma annihilation (which excludes electron-positron 
triplet state contribution and considers only singlet state annihilation). A^ is the number density of atoms or molecules 
in the medium; 5 is the Dirac <5-function; x and rj are the positron and the electron co-ordinates. 

Here we present a new formulation whereby the r.h.s of equation 1.2 is represented exactly by the onshell and 
half-offshell coupled-channel T-matrix scattering amplitudes, thus facilitating the evaluation of Z e jf without any 
explicit use of the scattering wave function. When a positron collides with a target, it can have direct as well as 
rearrangement scattering. So, we present our formulation in two different sections consisting of I) direct (elastic and 
inelastic) scattering and II) direct plus rearrangement (positronium formation) scattering. 



II. POSITRON ANNIHILATION CONSIDERING THE DIRECT SCATTERING CHANNELS 

In this section, we consider the direct scattering of positrons from atomic (<f> n ) targets where the total wave function 
is expanded as [9]: 

i>£ ( x > r i> r 2: ••• r Jv) = F n(x)^„(ri,r 2 , ...r N ) (2.1) 

n 

where F n 's are the expansion coefficients, representing the motion of the positron with momentum k; Vj is the 
co-ordinate of the j — th electron and x is that of the positron. The total Hamiltonian is partitioned as 

H = H° + V d (2.2) 

where H% is the unperturbed part of the total Hamiltonian in the direct scattering channel of the positron and the 
atom and V d is the positron-atom interaction potential. The unperturbed and the total Hamiltonians satisfy the 
following eigen-value equations: 

H° d \k<j> n ) = E\k<p n ) (2.3) 
(H° d +V d M+)=E\^) (2.4) 
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where E = k 2 /2m — £a is the total energy; Za is the binding energy of the initial target atom (<p n ); m and k are 
the reduced mass and the onshell momentum of the positron. With the eigcn-value equations 2.3 and 2.4 for the 
unperturbed and the total Hamiltonians, one can write the Lippmann-Schwinger equation for the scattering wave 
function \ip£) as [9,10]: 
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\W = + fl_ g o + j0 W) (2-5) 

Using the closure relation (2tt)- 3 $ d z k"\k" 4> n ,,)(k" 4> n „\ = 1, and using the T-matrix definition: Vd\tp^) = 
T\k(j> n ), in cqn(2.5), we arrive at the expression for the total scattering wave function in terms of the T-matrix 
elements: 

n" — 1 n 

The total scattering wave function can be evaluated from this equation [4]. However, we generally solve it for the 
scattering T-matrix amplitudes obtained by multiplying eqn.(2.6) by Va and projecting with (k'<f> n >\ and using the 
T-matrix definition: Vd\ip^) — T\k(f> n ): 

(k'<P„,\T\kcj> n ) = {k'(t> nl \V d \k(j> n ) 

1 f d 3 ] j / (k'<l>n'\Vd\k"cl ) ^)(k"^\T\k^ n ) 

E 1 Z7V/ Al\ \ ' ) 



(2n) 3 ^fJ E-E» + iQ 

Eqn.(2.7), in its one-dimensional partial-wave form (eqn.2.16), is exactly solved using the matrix inversion method 
[11]. Simultaneous equations are formed by replacing k' with various values of k" on which the radial integral over dk" 
is discretized. The solutions of the simultaneous equations give us the both-onshell ((k'(j> n '\T\k(j> n )) and half-offshcll 
({k"(j)n'\T\k(j) n )) T-matrix amplitudes for various values of k", where k, k' are on shell momenta and k" are the off- 
shell ones. While the solutions for the onshell T-matrix elements reflect the asymptotic behavior of the wave function 
and provide the physical cross sections, the half-offshell elements are usually thrown away. We understand that the 
latter might contain the short-range properties of the wave function and they together with the on-shcll elements can 
lead to an exact evaluation of the Z e ff. We multiply eqn.2.6 from left by J2f=i <^( x ~~ r j) = A, (say) and project it 
by (V'fc I to obtain: 



(^ + |A|^+) = (^+|A|fc^) 

+ J^f^j dk E-E^ + m (2 ' 



At this stage, to calculate Z e ff, which is equivalent to (^|A|^) (see eqn.1.2), we have two options: 1) using 
equation 2.6, substitute for (ip£\ in the r.h.s of (2.8) or 2) evaluate (ip£ \A\k<f) n ) separately and substitute in eqn(2.8). 
The first case leads to a complicated equation as follows: 

(V^|A|^> = (fc^|A|^„) 

+ f d3k ,,(k<l>n\T\k , 'cf )n// )(k ,, cj )n ,,\A\k ( l> n ) 



i^f^J E-E»-iQ 



4 4?xf*'* 



j3; „(A:^|A|fc"^)(fc'V„H7 1 |^„) 



(2irf^J E-E>> + i0 

^^fc^iTifc'/^^^'^HAI^C}^'^''^!^) 



(2.9) 

This equation, although can be solved numerically, needs an extra effort to evaluate the principal value part of the 
last term, which contains a product of two singularities arising out of the product of Greens functions. We, therefore, 
look for the evaluation of (k(j} n \A\ip^) by projecting eqn(2.6) with {k<f) n \A: 
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(k' ( t> n/ \A\i>+) = {k'4> n ,\&\k4> n ) 

+ J^2j] dk E-E» + iO (2 ' 10) 

We solve this equation exactly (which is a very straight forward numerical summation) and substitute the complex 
conjugate of (k'(/) n >\A\ip£) in cqn(2.8) to get Z e ff. However, like the T-matrix equation, we solve them in their 
one-dimensional partial wave form. To arrive at the corresponding partial wave equations for (2.10) and (2.8), we 
define the matrices D and Z as: 

A|V£ > = £>N„) (2.11) 
(^+\D = (k<j> n \Z (2.12) 

and rewrite eqn(2.10) and eqn(2.8) formally in terms of them: 

(k'(t>n>\D\k<t> n ) = (k'<t> n , \A\kci> n ) 

+ W?l?J E-E» + n (2 ' 13) 

and 

(k<p n \Z\k<p n ) = (k<f> n \D\k<t> n ) 



3; „ <^„|D|fc>„»}<fc"<MW«} 



where D is the complex conjugate of D. Using a partial wave decomposition of the form: 



(2.14) 



{k'fa \X\k<t> n ) =EEEEEE ( L ' l 'i M L ,m v \JM)Yi, Mv (k') 

J M L M L L' M l , 

(LI, M L m l \JM)Y LML (k)Xj(n'l'L'k', nlLk) (2.15) 

where, X = T, V, D, A, or Z; n, I are the principal and orbital quantum number of the target and L is the orbital 
quantum number of the moving positron in the initial state; primed quantities denote the same for the final state. 
With the above expansion, the scattering T-matrix equation and the above two equations for the D- and Z-matrices 
reduce to: 

Tj(T',k';T,k) = Vj(T',k';Tk) 



m " X- [ „ JllJ ,2 V,i{r',k'-T''X)Tj{r'' 1 k'' ] T,k) 

I^V ^fc» 2 + *o (2 - 16) 



Dj(T',k';T,k) = Aj(T',k';Tk) 
— V fdk 

4tt 3 2-f J " kl„ - k" 2 + iO 



+ s e / ^^ 2Aj(r ' fc :r' fc :?p (T ? v ' fc) 



Zj(T',k';T,k) = Dj(T',k';Tk) 

_/ !./._// l,ir\rr,l~ir l,". 



4- m " V f , 1? n h „2 Dj(T\k'-T"X)Tj{r"X;r,k) 
+ 4^ dkk k* T „-k»*+iO (2 - 18) 

T T 

where r = (nlL) and r' = (n'l'L 1 ); m" is the reduced mass of the projectile in the intermediate state (here, m" = m = 1 
in au). We suppress the suffix d from Vd for convenience. 

In terms of partial wave Z-matrices, Z e ff(k 2 ) comes out to be: 

Z eff (k 2 ) = V ——Zj(nlLk; nlLk) (2.19) 

z — ' 47T 

J 
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While eqn.(2.16) is generally used to study positron-atom scattering, eqs.(2.17) and (2.18) are particularly useful to 
evaluate Z e ff from the onshell and half-offshcll T-matrix outputs of eqn.2.16). We shall present a simple numerical 
account on e + -He scattering to verify the code and compare the numbers, but beforehand we present a general formula 
for Z e ff by inserting eqn.(2.17) into eqn.(2.18). This latter is of particular interest, as it will explicitly demonstrate 
how the Z e ff is dependent on the T-matrices. For this, we first rewrite eqs.(2.16), (2.17) and (2.18) in the following 
notations: 

Tk'k = Vk'k - iVk'k'Tk'k + Vk'k"G (k', k")T k "k (2.20) 
D k 'k = Afe/fe - iVk'k'Tk'k + Ak' P "G (k' ,p")Tp" k (2-21) 
Zk'k = D k 'k - iDk'k'Tkk + Dk' q "Go{k , q")T q n k (2.22) 

where summations over intermediate states are implied and off-shell momenta are represented by k", p" and q". In 
the above, we have used the following relation for the complex Greens function 

G+(fc 2 - k" 2 ) = \ = -inS(k 2 - k" 2 ) + ^— * (2.23) 

0V ' k 2 -k" 2 + i0 V ' k 2 -k" 2 y ' 

to expand it into real and imaginary parts; P represents principal value integration; Go represents the real (principal 
value) part of Gq . Inserting Dkk in eqn(2.22) we obtain an explicit relation for Z e ff(k 2 ) = Zkk'- 

Zkk = A/jfe + iAkk[Tkk — Tkk] + Tkp"GoA p »k + ^k q "GoT q »k 

+ &kk\Tkk\ 2 +T k p"G A p » q >iG T q » k (2.24) 
= A fefe + 2A fcfc Im[T fefe ] + A k k\Tkk\ 2 + T k p"G a A p n k + A kq "G T q nk 

+ Tkp"GoA p n q rrGoT q iik (2.25) 

where A pq corresponds to plane-wave value of Z e fj for the initial and final momenta p and q. lm[Tkk] and \Tkk\ 2 
are proportional to physical cross sections (representing the asymptotic behavior of the wave function) . Others are 
interference terms, linear and quadratic in the half-offshcll T-matrix elements, and expected to play a crucial role at 
low and intermediate energies. To understand their role, and to check the normalizations of equations 2.17 and 2.18, 
we provide a numerical test below. 



A. Numerical test to equations 2.17, 2.18 

To test the formulations of the equations 2.17, 2.18 (Equation 2.16 is well established in literature) we perform 
sample calculation on e+-He elastic scattering, considering only the static interaction [keeping summation over r" = 
He(ls 2 )], and evaluate the Z e ff using the resulting onshell and half-offshell T-matrix elements. We use atomic units 
throughout and use delta- function normalization for the plane wave. 

We note that equations {2.16,2.17,2.18} and {2.20,2.21,2.22} are equivalent. In figure 1, first we plot the dotted 
curve which is obtained considering only the plane wave parts (first term of the r.h.s) of eqs.(2.20), (2.21) and (2.22). 
This plane wave approximation gives a value of Z e ff — 2.0(= Z) as was expected and provides the normalization. 
Next we consider first two terms of the r.h.s of eqs.(2.20), (2.21) and (2.22) and plot the result as dashed curve. This 
approximation is equivalent of considering the plane wave and the onshell T-matrix contributions of equations (2.16), 
(2.17) and (2.18) leaving aside the half-offshell contribution (particular integral part of the Greens function). We 
obtain a lower value of Z e ff with increasing energies, signaling the manifestation of a repulsive potential at higher 
energies. Now, the solid curve is obtained with all the three terms of eqs.(2.20), (2.21), and (2.22). That is, considering 
both the onshell and the half-offshell contributions together with the plane wave. The static potential in a e+-Hc 
scattering is repulsive and consequently it lowers the value of Z e ff (and hence the annihilation rate). We compare the 
solid curve with the results of a Schwinger Multi-Channel (SMC) calculation on e + -He with the same physical content 
(considering static interaction only) [12]. Both the curves agree quite well. The marginal difference in the Z e ff value 
between the present and the SMC calculation is supposed to be acceptable since, the wave functions for Helium used 
in these two calculations are different (we use the Roothaan-Hatree-Fock five-term wave function of Clementi and 
Roetti [13] for He) and also the T-matrix and the SMC formulations are different. 

The results provided in figure- 1 are aimed at checking the new equations vis-a-vis other methodologies where we 
considered the static interaction. Formally, to arrive at a physically converged result for the Z e ff, full expansion 
basis indicated in equations 2.17 and 2.18 need to be employed like the case for the scattering T-matrix equation 
2.16. However, in practice it has been found that the positron-scattering cross sections do not converge easily unless 
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real and virtual effects of the rearrangement channel of positronium (Ps) formation is considered in the theoretical 
formulation. So, for a converged description of the scattering and annihilation we need to consider the Ps formation 
channel explicitly in the theoretical formulation. However, the above formulation is suitable for employing model 
polarization potentials alongside the static potential so as to arrive at a meaningful physically converged result 
without being confined in the ab initio framework. 



III. POSITRON ANNIHILATION CONSIDERING THE DIRECT AND THE PS-FORMATION 

CHANNELS. 

When the possibility of a real or virtual positronium formation is considered, through the capture of a target 
electron by the incident positron, the theoretical formulation for a single electron target differs from a multi electron 
target in the sense that for the latter case the Ps-target(ion) wave function need to be formally antisymmetrized. 
Here, wc discuss them in two different sections A and B. 



A. Single Electron Target 



For positron scattering from a single electron target, the total wave function (2.1) can be expanded (considering Ps 
formation) as [9,14]: 

^+(x, n) = ]T F„(x)^ n (n) + &(pi)x*(ti) (3-1) 

n v 

where p\ = (r\ + x)/2 and ti = ri — x. Q v and \ u represent the moving and the bound-state (z/-th) positronium 
atom. The total Hamiltonian is now partitioned as: 

H = H° d +V d = H° + V c (3.2) 

where H d , H® are the unperturbed Hamiltonians in the direct (d) and capture (c) channels satisfying the eigen- value 
equations 

H%\k<f> n ) = E n \k<f> n ) (3.3) 
H°\k Xv ) = E v \k Xv ) (3.4) 

and Vd and V c are the interaction potentials therein. E n = fc^/2 — £a and E v — kp s /A — £p s : £a and £p s are 
the binding energies of the initial target atom and the rearranged positronium atom; k x and kp s are the momenta 
of the positron and the positronium. In terms of the two-cluster channel- Greens-functions G d = (E — H®)" 1 and 
G° = (E — if°) _1 , we take the Lippmann-Schwinger integral equation for the wave function as [15]: 

|V+> = \kcj> n ) + G° d T d \kcj> n ) + G° c T c \k<p n ) (3.5) 

where, T d and T c arc defined as Vd\tp^) = T d \kcf) n ) (here T d = T of section-I) and V c \ip~£) = T c \k<f> n ). Using the 
following closure relations for the direct and the rearrangement channels, 

1 = pE/ *"|fc>n")(fc"<M (3-6) 

1 = pE / dk"\k" Xv „){k" X „„\ (3.7) 



we rewrite eqn.(3.5) as: 



. 3] „W'<l>n")(k"(pn»\T d \k<f> n ) 



(2 



W^J E-E'J + iO (d ' 8) 



G 



Here, ££' = k" z /2 - £ A and E'J = k" z /4 - £ Ps are the off-shell energies in the direct (d) and the capture (c) channels. 
We construct the coupled equations by 1) multiplying this equation with V d and projecting out with (£;'</>„< | and 2) 
multiplying this equation with V c and projecting out with {k'\ v i \ 



(2nf^J E-E>> + iO 

+ ^Z 7f J dk E-E'J + iO 



,3,// WXv' I Wtfn" \T d \k<t> n ) 



(k'XW \Wn) = {k' Xv , \V c \k4> n ) + -^3 ^ f d*k» E £ „ . ( 

1 \- f , 3 ,„ (k'x»> I Vc I fc"x„» ) {k" X v» \T C I ^ ) ._ in . 

+ j2^2L f J dk S-^' + zO (3 ' 10) 



where in eqn.(3.9), we use = -ff? + V c — (see eqn.3.2) and also use the eigen-value equations (3.3,3.4). 
Once the above coupled-equations are solved and we are equipped with the T-matrix amplitudes (p<j) n i \T d \k<f> n ) and 
(qXv' \T c \k(f) n ) for on-shell and off-shell values for the momenta p and q, we can get Z e ff in terms of them. To deduce 
Z e ff, in terms of T-matrices, we project equation (3.8) by (tp^\A and arrive at: 



(^ + |A|^ + ) = (^ + |A|fc^) + ^3^ J d 3 k" 



(i>+\A\k"cf> n „)(k"cf> n „\T d \k<p n ) 



E-E'^ + iO 

(2tt) 3 ^fJ U " E - £;' 



~ n (rt\A\k"x v ,,){k"xATc\k<j> n ) 



d 6 k n vrjo i-i- w is- ^ ^r-rn, ( g n) 



Now, {tp£\ from eqn(3.8) may be substituted in the r.h.s of eqn(3.11) to arrive at a direct expression for Z e ff. However, 
that will lead to a complicated equation like (2.9). We rather develop simpler equations to evaluate {tp^ \A\k(f> n ) and 
(V>fc |A|fcx„) and substitute them back in eqn(3.11). For this, we project eqn(3.8) from left by (k'(fi n >\A and (k'xv'\A 
and obtain: 



(fc>„,|A|^) = (fcVn'|A|fc0„} + -^3^ f d 3 k 

\ / /y-. II 



, l {k'ct> nl \A\k" ( j )nll ){k" ( j )nll \T d \k(jy n ) 



E-E" + zO 



1 \ " /',.;,.// " \- " \<-" I •■<■ ""„/ (■},-)■ 

■ WT^fJ '' " E-E'J + iO 



, 3 u , (k'K<\&\k"Xv>>){k"x,"\T c \kK} 
' + i0 

:\ u ll(k' Xv >\A\k" (p n » ) (fe"0n" \T d \k(j> n ) 



(k' Xvl \m) - (k' Xv ,\A\kK) +JL_£y>* E _ E ,, /() 



3u „(k'x„'\A\k"x v »)(k"x,"\T c \kci>n) 



(2tt)3 ^ J E-E'J + iO 



(3.13) 



The above two equations are very straight forward to solve as one need to carry only numerical integrations with 
known values of T d , T c and the calculated plane-wave matrix elements concerning A as inputs. We are not interested 
to repeat the calculations for T d and T c and rather hope that the existing T-matrix results [16] may be applied to 
calculate Z e ff. 



B. Many Electron Target 

For multi-electron targets the formulation is very much similar to that of section IIA, except few fundamental 
changes. Without repeating the whole thing, we thus mention here about the necessary changes. For the positron 
scattering from a multi-electron target, the capture channel need to be explicitly antisymmetrized and expressed as: 

^fe (x,ri,r 2 , ..rjv) = ^ F n (x)0„(ri, r N ) + Ai ^ ^(pi)x^(ti)^( r 2, -,r N ) (3.14) 

n v[i 
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where ip represents the residual target ion and Ai is the antisymmetrization operator, which antisymmctrizes electron 
1 with other target electrons. The initial target wave function <fi is supposed to be antisymmetrized implicitly. The 
total Hamiltonian is now partitioned as: H = H d + V d = #c(j) + ^c(j)' wnere ^c(j) an( ^ ^0) arc tnc unperturbed 
Hamiltonian and the Ps-target(ion) interaction potential in the capture channel of the positronium formation, with 
the j-th electron being attached to the positron. Accommodating the Pauli exclusion principle for the rearrangement 
channel, the Lippmann-Schwinger integral equation is now written as: 

IV£> = Mn) + G° d T d \k<f> n ) + A,G cU) T cU) \k^ n ) (3.15) 

where T d and T c are defined as V d \ip£) = T d \k<j) n ) (here T d = T of scction-I) and V c ^\^) = T c ^\kcf) n ). Using 
following closure relations for the direct and the rearrangement channels: 

1 = 7^)sJ2 / dk"\k"4> n „){k"4> n „\ (3.16) 

n" —1 

1 = (2^)3 EE / dk"\k'J X ^^)(k'J X ^^\ (3.17) 

and proceeding in a similar way, we represent the Lippmann-Schwinger equation 3.15 as: 

1 V- f j3lll \k"K")(k"^\T d \k^ n ) 



(2tt)3 E-E'^ + iQ 

1 \ f ^3 h n A 3 I k "Xv" <f^>> ) (k"Xv" ' \T c{ j) | k(j) n ) } 



The rest of the procedures are exactly similar to those described in section IIA and are not repeated here. 

In summary, we present a new ab initio methodology to calculate Z e ff from physical (onshell) and virtual (half- 
offshell) scattering T-matrix amplitudes, without any use of the scattering wave function. The formulation presented 
here is for positron annihilation in atoms, but it could be universally applied to other annihilation studies as long 
as the dynamics of the interacting particles (or clusters) can be described by the well-known Lippmann-Schwinger 
type equation. The methodology is exact and can act as an useful tool for the annihilation studies as most of 
the scattering theories (T-matrix, K-matrix, S-matrix) yield directly the scattering amplitudes. Performing elastic 
scattering and employing clastic-channel T-matrix amplitudes (on- and off-shell) we reproduce the corresponding Z e f f 
result and demonstrate the utility of the methodology. A similar T-matrix formulation for the pick-off annihilation 
of ortho-positronium is under way and will be published soon. 

The work has been carried out under financial support from FAPESP, Brazil through project number 99/06844-7. 
I gratefully acknowledge various discussions with Dr. T. Frederico, Dr. J. S. E. Germano of our department and Dr. 
M. A. P. Lima and Mr. M. Varella of UNICAMP, SP, Brazil. 
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Figure Caption: 

Figure 1. Theoretical values of Z e ff in various approximations as a function of positron energy for the target of 
atomic helium. 
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